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Abstract 

The 7] production in the {n,n') bottomonium transitions T(n) — > 
T{n')r], is studied in the method used before for dipion heavy quarko- 
nia transitions. The widths r^(n, n') are calculated without fitting 
parameters for n = 2,3,4, 5, n' = 1. Resulting r^(4, 1) is found to be 
large in agreement with recent data. Multipole expansion method is 
shown to be inadequate for large size systems considered. 

1 Introduction 

The rj and 7r° production in heavy quarkonia transitions is attracting atten- 
tion of experimentahsts for a long time pj. The first result refers to the 
ip{2S) J /ip{lS)r] process (to be denoted as ip{2, l)r] in what follows, sim- 
ilarly for T) with ^ = (3.09 ± 0.08)% [1], Ttot = 337 ± 13 keV. 

For the T(2, 1)?7 and T(3, l)ri transitions only upper limits B < 2 ■ 10~^ 
and B < 2.2 ■ 10"'^ were obtained in [2]| and [3j correspondingly and prelim- 
inary results appeared recently in |5, i?(T(2, l)ri) = (2.5 ± 0.7 ± 0.5)10""^, 
5(T(2, l)7rO) < 2.1 • 10-^(90% c.L) and E(T(3, l)r/) < 2.9 • 10"^ in [S]. 

On theoretical side the dominant approach for both dipion and single rj 
and TT production is the Multipole Expansion Method (MEM) (see [6l [7j and 
refs. therein), where it is asumed that heavy quarks emit two gluons and the 
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latter are converted into meson(s) by a not clarified mechanism. An essential 
requirement for this mechanism in QED is the smallness of the source size tq 
as compared to the wavelength, so that rok -C 1. 

In reality both for charmonium and bottomonium transitions rok > 1, 
but it is not this parameter which invalidates MEM for heavy quarkonia. It 
appears, that in QCD there is another important length parameter, the QCD 
vacuum correlation length A, which makes it impossible to emit freely gluons 
at points separated by distance r, r > A. 

The value of A was found on the lattice and analytically, A < 0.2 fm 
[8]. Since r.m.s. radii of all excited cc,bb states are larger than 0.5 fnflall 
vacuum gluons there are correlated forming the QCD string and emission 
of additional gluons (if any) implies formation of heavy hybrids. All this is 
considerd in detail in Field Correlator Method (FCM) [9j. 

One can make an independent check of MEM in application to the bot- 
tomonium level calculation. Here MEM yields nonperturbative correction to 
the levels expressed via gluonic condensate [10]. Comparison to the exper- 
imental data shows (see [TT] and Table 1 below) that for all level splittings 
except (25* — IS") in bottomonium, MEM prediction is more than 50% off, 
while for charmonium MEM does not work at all. Thus one concludes that 
only at distances below or equal 0.2 fm, MEM can give reasonable results, 
while for all states of charmonium and all excited states of bottomonium 
(where sizes are much larger than vacuum correlation length A) the applica- 
tion of MEM is unjustifiable. 



Table 1. Predicted splittings, between spin- averaged levels, and experi- 
ment. 



splitting 


MEM[U] 


FCM[12] 


exp.[l] 


2S - IS {bb) 


479 


557 


558 MeV 


2S - IP (bb) 


181 


122 


123 MeV 


3S - 2S (bb) 


4 570 


348 


332 MeV 


2S - IS (cc) 


9 733 


610 


606 MeV 



"'^In reality the r.m.s. radius tq of heavy quarkonia is not small, e.g. for T(riS) it is 
equal to 0.2 fm, 0.5 fm, 0.7 fm, 0.9 fm, for n = 1,2, 3, 4, 5, and for charmonium this radius 
is even larger: 0.4 fm and 0.8 fm for n = 1 and 2 respectively. 



2 



A similar failure of MEM is found in applications to dipion bottomonium 
transitions, where using MEM one can fit dipion spectra in T(2, 1)7111, but 
not in T(3, l)7r7r and T(4,2)7r7r |6l [7j. In contrast to that, FCM as will be 
discussed below explains both spectra and cos 6* dependence for all dipion 
transitions in universal approach with two fixed parameters. 

In FCM large distances are under control and not single gluons but com- 
bined effect of all gluons in the string defines the dynamics. 

In particular, single eta emission in heavy quarkonia proceeds via string 
breaking due to qq pair creation with simultaneous emission of vr or t]. The 
flavor SU(3) violation in 77 production then resides in difference of threshold 
positions and wave functions for BB and BgBs {DD and DgDs) intermediate 
states. 

As it is clear, the dynamics of FCM for rj emission does not depend 
strongly on heavy quark mass, and only sizes of initial and final heavy quarko- 
nia states and intermediate heavy-light mesons enter in the form of overlap 
matrix elements. 

In contrast to that, MEM predicts a strong dependence on the heavy 
quark mass, r^(?7,,n') = O(^). In addition in [7] a strong suppression 

of the ratio T^/T.,^t^ with the growth of the energy release AM = M(n) — 
M{n'), F^/F^^ ~ (AM)7 predicted for higher excited states of quarkonia 
and bottomonia, which does not agree with experiment (see below). 
Using models based on MEM in [6] small ratios of widths 

F(T(2,l)r/) , , F(T(3,1)?7) , 

^) , ; ' i '[ ^ 2.5 ■ 10"^ and ^) ! '[ = 1.3 ■ 10"^ 1 

F(^(2,l)r7) r(^(2,l)r/) 

have been predicted, with the model property that the bottomonium yields 
of r] would be smaller than those of charmonium; specifically in the method 
of [6], the width is proportional to p^/ttlq, so that for T{4,1)t] the ratio 

r(T(4,l)r;) _ o o in-S 

However recently [13] new BaBar data have been published on T(4, l)r) 
with the branching ratio 

5(T(4, l)ri) = (1.96 ± 0.06 ± 0.09)10^^ (2) 

and [13] 

F(T(4,1)77) 

K K , J u ^ 2.41 lb 0.40 ±0.12. (3) 



F(T(4, l)7r+7r- 
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This latter result is very large, indeed the corresponding ratio r(vi(2'ijr;) 
is ~ 0.4 and theoretical estimates ([H) from [6j for a similar ratio yields 
3.3 ■ 10~^.Thus, the experimental ratio is very large as compared to MEM 
predictions [6l[7j. All this suggests that another mechanism can be at work 
in single r] production and below we exploit the approach based on the Field 
Correlator Method (FCM) recently applied to T(n, n')7r7r transitions with 
n < 3 in [HIlS], n < 4 in [16j and n = 5 in [HI IE]. 

In this paper we confront MEM and FCM and show that recent exper- 
imental data on single r] production in T(45') — T(IS') transition give a 
strong support to the FCM result and cannot be explained in the framework 
of MEM. 

The method essentially expoits the mechanism of Internal Loop Radiation 
(ILR) with light quark loop inside heavy quarkonium and has two fundamen- 
tal parameters - mass vertices in chiral light quark pair qq creation Mbr ~ /tt 
and pair creation vertex without pseudoscalars, ^ 2u!, where uj^uJs) is 
the average energy of the light (strange) quark in the B{Bs) meson. Those 
are calculated with relativistic Hamiltonian [12] and considered as fixed for 
all types of transitions u = 0.587 GeV, Us = 0.639 GeV (see Appendix 1 of 
[I4j for details). 

Any process of heavy quarkonium transition with emission of any number 
of Nambu-Goldstone (NG) mesons is considered in ILR as proceeding via 
intermediate states of BB,BB* + c.c.,BsBs etc. (or equivalently DD etc.) 
with NG mesons emitted at vertices. 

For one rj or 7r° emission one has diagrams shown in Fig.l, where dashed 
line is for the NG meson. As shown in [T5| [16] . based on the chiral 
Lagrangian derived in [191, meson emission vertex has the structure 







<fa\ 





71 



K+ \ 



[/ = exp 



•a 



n 



(5) 



4 




Fig.l Single eta production (dashed line) from T{n)BB* vertex (a), and BB*T{n') 
vertex (b). 

The lines (1,2,3) in the tj matrix ^ refer to u,d,s quarks and hence to the 
channels B^B^ , B^B^, B'^B^ (and to the corresponding channels with B* instead 
of -B). Therefore the emission of a single r] in heavy quarkonia transitions requires 
the flavour SU{3) violation and resides in our approach in the difference of channel 
contribution BB* and BsB*, while the vr*^ emission is due the difference of B^B^* 
and B^B^* channels (with B ^ D for charmonia). 

The paper is devoted to the explicit calculation of single rj emission widths in 
bottomonium T(n, l)r/ transitions with n = 2,3,4,5. Since theory has no fitting 
parameters (the only ones, and Mbr are fixed by dipion transitions) our pre- 
dictions depend only on the overlap matrix elements, containing wave functions 
of T(nS'), B, Bg, B* , B*. The latter have been computed previously in relativistic 
Hamiltonian technic in [12] and used extensively in dipion transitions in |161 I17 1 [T8]. 

The paper is organized as follows. In section 2 general expressions for process 
amplitudes are given; in section 3 results of calculations are presented and discussed 
and a short summary and prospectives are given. 

2 General formalism 

The process of single NG boson emission in bottomonium transition is described by 
two diagrams depicted in Fig.l, (a) and (b) which can be written according to the 
general formalism of FCM |14| , I16| , \T7\ as (we consider rj emission), see Appendix 
for more detail, 

M = + A^f ); = A^gij; " ^bb* > ^ = 1> 2 (6) 
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For the diagram of Fig. 1(a) the amphtude for intermediate BB* or B^B* state 
can be written as 

i')(p,k)J„,(p) d^p 



(7) 



^-^(p) (27r)3' 

while A^r; , corresponding to the diagram of Fig. 1(b), has the same form, but 
without NG boson energy in the denominator of ([7]). The overlap integrals of 
T{nS) and BB* wave functions with emission of r] with momentum k are denoted 

by Jn\p,^), the corresponding integrals without rj emission are given by Jn'(p)- 
Finally we define all quantities in the denominator of ([7]); in A^^^* the de- 
nominator is 

E-E{p) = M{T{nS))-{u;,+MB+M*s+^+^^-^) ^ -AM*-a;,-S(p, k). 



2Mb 



2M% 



(8) 

For Mil one omits ujri and k in ([8]). Finally after taking Dirac traces in amplitudes 
and accounting for the p-wave of emitted r/ one can represent the matrix element 
M.r) as follows: 



1 



UJ 



3 s 



1 



£(1) 



(9) 



Indices i'i in Ci^u in Q refer to the T{n'S) and T{nS) polarizations respectively. 



(10) 



with £^'^\C^s^ being integrals of overlap factors J„(p, k) J„/(p). 

The width of the T(n,n')r] decay is obtained from |A^p averaging over vector 
polarizations as 



(1) 



+ 



(2) 



£(2)' 



(11) 



where the phase space factor d<I> = p^27r5(M(T(n))— M(T(n'))— w^— 2]g^Y(;^' 
Introducing the average a) = \{uJs + 1^), one can rewrite (fTT|) as 



2w 



fc3 



2^ 



UJ 



(12) 
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0.063. 



with C, = ^^^ 

One can see from the general structure of F^, that the main effect comes from 



the difference 



10.882 - 1.139| PS 0.257, and from the difference of 



3 Results and discussion 

We consider here the single rj emission in bottomonium transitions T(n, l)rj with 
n = 2, 3, 4, 5. The corresponding values of AM*, AM*, w,,, k are given in the Table 
2. 

Table 2. Mass parameters of T{n,n')r] transitions (all in GeV, k in GeV/c). 



(n, n') 


2,1 


3,1 


4,1 


5,1 


AM* 


0.582 


0.25 


0.026 


-0.26 


AM* 


0.757 


0.425 


0.20 


-0.08 


ujn 


0.562 


0.87 


1.075 


1.325 


k 


0.115 


0.674 


0.923 


1.20 



The resulting values of r^(n, n') have been computed as in (I12p with to = 0.587 
GeV and lJs = 0.639 GeV, calculated earlier in [12], see Table 4 of [Hj, and with 
wave functions fitted to the realistic wave functions in [16], (set I), while in set II 
parameters of the T{nS) wave function were changed by 10-15%. 

Results of calculations are given in Table 3, where we have put M^r ~ fn and 

(^) ~ ^ to explain the decay Tbb{4:S) = T^xp (see [16] for details). 



Table 3. Values of r.^(n, n') (in keV) calculated using Eq. ^ with (^|^ 
^ vs experimental data r°^P(n, n') (in keV). 



(n, n') 


(2,1) 


(3,1) 


(4,1) 


(5,1) 


r^, set I 


2.50-10^^ 


1.45 


0.9 


3.5 


r^, set II 


1.4-10"^ 


3.6-10"^ 


1.9 


1.3 


n 


(0.8±0.3) • 10"^ 

[4] 


< 5.8 • 10"^ 
[5] 


4.0 ±0.6 

[13] 
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Looking at the Table 3, one can see, that there is an order of magnitude 
agreement with experiment for the set II. Indeed, the factor (^y^) ('^) 
estimated from T(n, n')7r7r transitions studied in [I4j-[18] to be roughly in the range 
[^,2]. On theoretical side our formulas ([7|)- pT]) automatically produce the width 
r^(n, n') of the order of 0(1 keV), for all (n, 1) transitions except for (2,1), where a 
small phase space factor k'^ gives two orders of magnitude suppression of r^(2, 1), 
and for r^(3, 1), which is highly sensitive to the form of wave function (cf. set I 
and set II). For the r^(5, 1) one obtains a 7 keV value, which is however small as 
compared with the r7r7r(5, 1), the latter being 0(1 MeV). For r^(4, 1) and r7r7r(4, 1) 

from [16] the calculated ratio is i?r,/7r7r = r'^^'lti) ~ ^ ('^) (l^) ^ which 
roughly agrees with experimental value R'^^^ = 2.41 it 0.40 ± 0.12 |13| . 

To compare our results with MEM predictions and experiment, we list in Table 
4 the ratios r(V'(2,'i)^) ^ ^ ^ ~ 2, 3, 4, 5 (the MEM numbers for n = 4, 5 are 
obtained from n = 3 using scaling X = 0{p^) |6J. 

Table 4. Values of the ratio rSja.'ih) = X for n = 2,3,4,5 from MEM [6], 
experiment and the present paper (for the latter we take rexp('0(2, 1)^) in the 
denominator and two sets from Table 3): 



n 


2 


3 


4 


5 


MEM 


0.25 • 10"^ 


0.13 • 10"^ 


0.33-10-^ 


0.73 • 10"^ 


exper. 


(0.8±0.2) • 10"^ 


< 0.58 • 10-^ 


0.40 ± 0.06 






[4] 


[5] 


[13] 




present paper 


(1.4^2.5) • 10~^ 


(0.36- 10^^^0.15) 


0.9^ 1.9 


0.13^0.35 



To check stability of our results, we have used for the wave function of Bg 
the realistic wave function different from that of i3. As a result one obtains for 
, ^''^"'^^ the values (2.74 • 10"^; 1.13; 0.44; 7.3) keV for n = 2,3,4,5 respec- 

tively. To understand why r.^(.l) is so sensitive to the form of the wave function (cf. 
results for the sets I and II), we have varied the shape of wave functions of T{nS) 
with n = 2, 3, 4, 5 (explicitly the parameter f3i) in the range ^ = ^^pt"^^ = 0.8 H- 1.2 

of the optimal value, reproducing the realistic wave function. We have obtained 
the dependence of r^(n, 1) on ^ for n = 2, 3, 4, 5. As a result r^(2, 1) is rather sta- 
ble, whereas r^(3, 1) has minimum (almost zero) for ^ = 1.17. This fact explains 
results of sets I and II; while set I produces normal (and nearly maximal) values 
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for r^(n, 1) of the order of 1 keV for n > 2, in case of set II the parameter /3i of 
T(3S') roughly corresponds to the minimum of F^. In this way doing comparison 
with experiment one obtains possibly an instrument for a precision study of wave 
functions of excited bottomonium states. 

Summarizing, we have calculated the single rj production width T^{n,n') for 
T(n, l)r] transitions with n = 2,3, 4, 5. We have found that r^(4, 1) are of the order 
of and larger than Tt^tj{4, 1). This fact is in agreement with the latest measurements 
in p] of r^''P(4, 1). We have shown that r^(n, 1), n = 4, 5 is large (~ 0(1 kev)) for 
typical (realistic) parameters of T(nS') wave functions, but can occasionally drop 
near zero for slightly varied form of wave function, as it happens for n = 3. This 
high sensitivity is connected to the oscillating character of excited bottomonium 
wave function^. Our calculations do not contain fitting parameters; the only 
two parameters M^j. , are fixed by previous comparison with dipion data. One 
should stress that rj production in bottomonium is not suppressed in our approach 
as compared to rj production in charmonium transitions. This is in contrast with 
the results of method of [6, 7j. We have given arguments why the dipion transitions 
in high excited states of heavy quarkonia as well as single rj and tt emission cannot 
be reliably calculated within the MEM method of [H [7], widely used now. As it 
is seen in Tables 3 and 4 the sequence of experimental data [1],[5], [E] contradict 
predictions of MEM. Recently a new calculation was done of T^{n,n') in |20j 
where also BB* etc. intermediate states were taken into account as well as in our 
approach. The authors however did not use wave functions of hadrons involved, 
but rather exploited fitted coupling constants and formfactors, and specific form 
of matrix elements, which makes it difficult to compare with our method directly. 

The authors are grateful to S.I.Eidelman for constant support and helpful 
comments. 
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4961.2008.2 is gratefully acknowledged. 

References 

[1] W.-M. Yao et al. Particle Data Group, J. Phys., G 32(2006) 1. 
[2] V.Fonseca et al., (CUSB Collaboration), Nucl. Phys., B 242(1984) 31. 
[3] I.C.Brock et al., Phys. Rev., D 43 (1991) 1448. 
[4] D.Kreinick et al., (CLEO Collaboration), arXiV: 0710.5929. 
■^A similar phenomenon of cancellation in matrix elements was mentioned in [6] 



9 



[5] Q.He et al, (CLEO collaboration'! . larXiv:0806.3027l 

[6] T.-M.Yan, Phys. Rev., D 22 (1980) 1632; 

Y.-P.Kuang, Front. Phys. China, 1(2006) 19 , arXiv: | hep-ph/060T044l 

[7] M.B.Voloshin, arXiv: 0711.4556; M.B.Voloshin, V.I.Zakharov, Phys. Rev. 
Lett. 45 (1980) 688. 

[8] M.Campostrini, A.Di Giacomo and G.Mussardo, Z.Phys., C 25 
(1984) 173; 

A.Di Giacomo, and S.Olejnik, Z. Phys., C 34 (1986) 577; G.Bah, N. Bram- 
billa, A.Vairo, Phys., Lett. B 421 (1998) 265. 

[9] for a review see A.Di Giacomo, H.G.Dosch, V.I.Shevchenko, Yu.A.Simonov, 
Phys. Rept. 372 (2002) 319. 

[10] H. Leutwyler, Phys. Lett., 98B (1981) 447; 

M. B. Voloshin, Nucl. Phys., B 154 (1979) 365; Sov. J. Nucl. Phys., 36 (1982) 
143. 

[11] Yu.A.Simonov, S.Titard, F.J.Yandurain, Phys. Lett. B 354 (1995) 435. 

[12] A.M.Badalian, A.I.Veselov, B.L.G.Bakker, J.Phys. G 31(2005) 417 
Phys. Atom Nucl. 70 (2007) 1764 . 

[13] N.Arnaud (for BaBar Collaboration), arXiv: 0805.2591. 

[14] Yu.A.Simonov, Phys. Atom. Nucl. 71 (2008) 1048, arXiv: 0711.3626. 

[15] Yu.A.Simonov, JETP Lett., 87 (2008) 147, arXiv: 0712.2197. 

[16] Yu.A.Simonov, A.I.Veselov, arXiv: 0804.4635; Phys. Rev. D (in press). 

[17] Yu.A.Simonov, A.I.Veselov, Phys. Lett. B 671 (2009) 55, larXiv:0805.4499l 

[18] Yu.A.Simonov, A.I.Veselov, JETP Lett. 88 (2008) 79, arXiv: 0805.4518. 

[19] Yu.A.Simonov, Phys.Rev. D 65 (2002) 0940018. 

[20] Ce Meng, K.-T.Chao, Phys. Rev. D 78 (2008) 074001. 

Appendix 1 

Matrix element of single r] emission 
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According to the general theory in |14[ 116] . the matrix elements \ Mc^^ 

for T(n, n')r/ corresponding to diagrams of Fig.l, (a) and (b) respectively, can be 
written as 

Here 77-2, ns are channels of intermediate state, with e.g. n2 = B, B* , Bg, B*, we 
omit indices n2 , and write 

Jn(p,k)= yy(^^23(^^n(cp- |+q)^„2(q)V^n3(q-k) (A.2) 

where ^rni^ni are momentum space wave functions of T{nS) and B{B*) mesons 
respectively. 

The vertex factor y^gs is calculated in the same way as in |14j . namely from 
the Dirac trace of the projection operators for the decay process in our case this is 

T(nS) — > BB*rj. Identifying the creation operators as 'i/';,7jV'b) V'ftTsV'n! V'bTfcV'nj ^ = 
u, d, sa and extracting vertex of r/ creation from the Lagrangian AC = — J •ipntJMffrtpnd'^x 

I 1 

1 I , one has for the decay pro- 

-2 

cess (cf. Appendix 1 of 

G(T BB*r]) = tr[-fiSb{u, w)j5Sn{w, x)j5Sn{x, w)-ikSi{w ,v)] (A. 3) 

As shown in [14J, appendix 1 and 2, the quark Green's functions can be 



which gives iM6rj^2^iMi^ with A = ^ 

fn^2ulr, V3 



split into two factors S{x,y) = A^G{x,y), with the projection operators A^ 



^^2ul^^^ ' ^ = n and the scalar part G{x,y), where spins are present only 
in spin-depdndent interaction and treated as corrections. Here is the average 
energy of quark in given meson. Hence one is brought to the spin factor Z. 

Z = tr(7,A+75A;75A+7fcA,-) (A.4) 

which is equal to 

z = ^^^i^ii-^' - - pIpI +pM)- (A.5) 

Here i},uj are average energies of b and n quark in B or B* . One can identify 
the momenta of B and B* as Pi = p and P2 = — p — k, then q in ()A.2p can be 
expressed as 
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and Z is (we put = 0) 

= ^(-kq(5ifc + /ci^fc + ^fc^i)- (A. 7) 

It is important, that we are looking for the P-wave of emitted r/, and hence for P 
wave of relative BB* motion, hence the integral ()A.2p should yield the term pk. 
This indeed happens, when one approximates ^n,tpn as series of oscillator wave 
functions and ()A.2p has the form 

Up,k)=yl,,e ^ W/„(p). (A.8) 

In the process of dq integration in ()A.2p one changes the integration variable 
qi ^ q[ — Un + 0{ki) with = (32, are oscillator parameters, found by c/ii^ 
procedure. 

Thus result of d'^q integration yields 

XL 

yi23 = ^(-kpf^ifc - hPk + ^fePi)- (A.9) 
In an analogous way one obtains for Jn'(p) in (|A.ip the form 

-Pi 

^'(P)=yi'23e ^'^4'(p) (A.IO) 
and yjj?23 obtained from the Dirac trace for the process BB* T(n'5'), 

Z2{BB*) = ^ei,ki{-2qi + -^Pi) (A.ll) 

and the result of integration over dPq yields in (|A.10p 

VuL = -^^'M-. (A.12) 

UJ 

Here i' is the polarization of T{n'S) (represented by 'ipbli''^b) and k as in (|A.9p 
is the polarization of B* . Averaging over angles of p one obtains 

,2 



and finally one writes as in 



fny/2LJr,V^ • 3 

For n {n,n') one can use time inversion and interchange indices and 
change sign of k, obtaining in this way Eqs. ([9]) and (llOp of the main text. For 
the intermediate state of B*B* the summation over polarizations of B* yields a 
net zero result, therefore we are left with only {BB* + B*B) intermediate state. 
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